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ABSTRACT 


Self-focussing  and  self-trapping  of  light  rays  are  studied  in  a classical 
two-dimensional  nonlinear  model  both  in  the  geometrical  optics  approximation 
and  when  retaining  higher  dispersion  terms.  The  paraxial  ray  approximation  of 
Akhmanov  (1)  that  leads  to  the  nonlinear  Schrodinger  equation  is  specifically 
avoided.  It  is  found  that  in  the  geometrical  optics  approximation,  focussing 
must  eventually  occur  once  convergence  of  rays  begins.  When  both  diffraction 
and  nonlinearity  are  present,  however,  focussing  may  or  may  not  occur. 
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SIGNIFICANCE  AND  EXPLANATION 


In  the  classical  theory  of  light,  rays  are  straight  except  for  the  smearing 
effect  of  diffraction  that  occurs  when  there  is  a sudden  change  in  intensity 
such  as  would  occur  when  light  is  partially  blocked  by  an  obstacle.  Then  light 
diffuses  to  some  extent  into  the  shadow  zone  so  that  the  boundary  between  light 
cuid  dark  is  not  sharp.  An  additional  effect  that  was  discovered  when  lasers 
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were  developed  is  that  in  some  media,  when  the  light  intensity  is  sufficiently 
great,  the  behavior  of  light  depends  on  its  brightness.  This  effect  is  called 
nonlinearity  because  the  governing  equations  beccme  nonlinear  and  one  of  the 
consequences  is  that  beams  of  light  can  be  caused  to  focus,  giving  such  great 
intensities  as  to  damage  the  crystal  medium.  Short  of  this  the  rays  can  become 
trapped  so  that  they  remain  in  a tight  bundle  and  never  spread.  This  paper 
studies  the  process  of  self -focussing  emd  self -trapping  on  a simple  model  and, 
in  particular,  refutes  the  assertion  of  Akhmanov  (1)  that,  for  this  model,  once 
rays  begin  to  converge,  a focus  must  result.  It  turns  out  that  the  smearing 
effect  of  diffraction  can  prevent  focussing  although  it  need  not  do  so. 
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The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  authors  of  this  report. 


SELF-FCCUSSED  AND  SELF-TRAPPED  BEAMS  IN  NONLINEAR  OPTICS 


Mohamed  N.  Y.  Anwar  and  Robert  D.  Small 


Introduction 


Light  rays  in  certain  nonlinear  media  converge  solely  due  to  the  presence 
of  the  beams  themselves.  The  property  responsible  for  self -focussing  and  self- 
trapping  of  beams  is  that  the  phase  velocity  of  the  waves  decreases  with  in- 
creasing amplitude.  A disturbance  that  has  a distribution  of  amplitude  along 
an  initially  flat  wave  front  will  be  retarded  at  its  most  intense  point,  as 
shown  in  Figure  la,  and  rays  will  bend  inwards.  The  convergence  of  rays  in- 
creases the  amplitude,  enhancing  the  nonlinear  effect,  resulting  in  greater 
convergence  until  a singularity  occurs.  Rays  bend  not  only  due  to  nonlinearity 
but  also  due  to  diffraction  which  resists  sudden  changes  in  amplitude  causing 
rays  to  bend  outweurds  to  disperse  the  energy.  When  focussing  begins  diffraction 
increases  to  oppose  it.  Diffraction  can,  in  fact,  restore  a focussing  beam  to 
its  original  state  as  in  Figure  lb.  In  this  paper  we  demonstrate  that  the  effect 
of  nonlinearity  alone  causes  focussing  but  that  the  combined  effect  of  nonlinearity 
and  diffraction  Ceui  lead  either  to  a focus  or  not,  depending  on  the  initial  con- 
figuration of  the  rays. 

The  analysis  for  typical  nonlinear  media  is  very  difficult  since  the  media 
are,  in  general,  anisotropic  and  also  require  a quantum  mechanical  description. 

The  present  work  studies  the  self-focussing  effect  of  a classical  isotropic  model 
in  two  dimensions.  This  has  been  attempted  before,  for  example  (1)  and  (2)  but, 
for  the  most  part,  previous  work  in  both  geometrical  and  dispersive  approximations 
has  depended  upon  a further  approximation  where  rays  must  be  nearly  parallel  to 
the  axis  of  the  beam  and  in  the  dispersive  case  the  governing  equation  is  the  non- 
linear Schrodinger  equation  (1) . In  the  region  of  any  focus,  however,  the  ray 
angles  necessarily  become  large  and  violate  that  approximation.  We  shall  show  by 
explicit  solutions  that  in  the  absence  of  diffraction,  nonlinearity  causes 
focussing  and  that  when  diffraction  is  present  it  can  prevent  focussing.  This  is 


in  contradiction  to  Atchnanov  et  al.  (I)  who  claim  that  diffraction  cannot  prevent 
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These  scalars  will  be  slowly  varying  as  well.  The  result  of  substituting  [3] 

into  [1]  and  [2]  and  equating  coefficients  of  cos  6 and  sin  8 , retaining  one 

order  of  smallness  in  a , is  the  set  of  equations 
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We  have  used  k = |k|  and  subscripts  to  indicate  partial  differentiation.  In 
addition  the  consistency  relations  from  [4]  are 
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We  now  consider  only  time  independent  envelopes.  Neglecting  time  derivatives 

in  (6]  we  see  that  u is  constemt  emd  neglecting  them  in  [5]  we  can  take  b = 0 
to  satisfy  the  last  equation.  The  third  equation  of  [5]  becomes  algebraic  and  is 
used  to  eliminate  b from  the  first  two  of  (5)  to  one  order  in  a . These  two 


remaining  equations  and  the  survivor  of  [6]  beoome 
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Equations  of  this  structure  have  been  seen  before.  Cumberbatch  (2)  and 
Shvartsburn  (3)  derived  them  for  a non-inertial  medium  and  hence  their  interpreta- 
tion of  [7]  differs  from  ours.  In  addition  Shvartsburg  used  an  exponential 
rather  than  cubic  nonlinearity.  We  emphasize  that  the  only  approximation  so  far 
is  the  neglect  of  higher  powers  of  a and  while  a,  and  must  be  slowly 

veurying  we  have  not  yet  neglected  any  of  their  derivatives. 

We  can  make  some  observations  immediately  by  solving  [7]  for  the  phase 


velocity  m/k  , obtaining 
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Neglecting  the  second  derivatives  and  the  term  in  a we  have  the  phase  velocity 

of  a uniform  plane  wave  in  a linear  medium.  If  = 0 the  medium  becomes 

P 

uncoupled  from  the  electromagnetic  wave  and  the  wave  passes  at  speed  c . With 
coupling  due  to  finite  its  speed  is  less  than  c and  depends  on  the  wave 

frequency  u . This  is  called  a dispersive  medium  because  the  group  velocity 
3(d/3k  depends  on  u so  that  wave  packets  of  different  frequencies  separate. 

When  (1)  = ojjj  the  wave  resonates  with  the  medium  and  cannot  pass  through  it. 

/~2  ^ 

Also  there  is  a frequency  band  U-  < u < /id.  + u where  the  wave  reflects  back- 

0 — — 0 p 

wards  from  the  medium,  the  amplitude  falling  exponentially  as  the  wave  penetrates. 

When  the  second  derivatives  of  amplitude  become  large  relative  to  the  ampli- 
tude there  is  an  additional  influence  on  the  phase  velocity.  These  derivatives 
are  called  higher  dispersion  terms  since  they  modify  the  previously  discussed 
dispersive  effects  emd  they  are  also  responsible  for  diffraction.  Near  the  peak 
of  a local  maximum  in  amplitude  we  have  a^  < 0.  This  negative  term  increases 
the  phase  velocity  relative  to  other  portions  of  the  wave  front  giving  a defocus- 
sing effect.  Finally  the  term  containing  a gives  the  nonlinear  influence  of 
the  medium  on  the  wave.  For  a > 0 this  term  causes  the  phase  velocity  to  fall 
as  amplitude  rises  so  that  a local  maximum  in  amplitude  tends  to  focus.  When 
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diffraction  and  nonlinearity  are  in  perfect  balance  we  obtain  the  self-trapped 
beam  of  Chiao  et  al  (4)  where  rays  are  straight. 

To  prepare  for  solving  the  equations  we  introduce  the  scaled  variables 


X*  = Fx  , y'  » Fy  , ic'  » 


Then,  if  the  primes  are  dropped,  181  and  [9]  do  not  change  form  while  [7] 
becomes 

[10]  * _ + « - - Da  + 2a^  - 0 . 


2.  Ray  Co-ordinates  and  Geaoetrical  Optics 


We  wish  to  investigate  solutions  to  the  systeo  of  equations  18] , [9]  , [101  • 
In  this  section  we  consider  as  negligible  the  terms  a^  and  euid  call  this 

the  geometrical  optics  approximation.  Then  in  Section  3 we  treat  these  higher 
dispersion  terms  equally  with  the  other  terms.  Equations  [B]  and  [9]  suggest, 
by  their  similarity  to  the  equations  of  conservation  of  mass  and  null  vorticity 
in  fluid  flow  problems,  the  introduction  of  the  potentials 
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since  those  equations  would  be  satisfied  identically.  Instead  of  eliminating  a 

and  K in  favour  of  ( and  n , however,  we  use  C and  n as  new  independent 

variables.  Then  the  vectors  along  co-ordinate  lines, 

2 

VC  • (Kj^rKj)  f Vn  = a (-|C2,Kj^) 

are  perpendicular  so  that  the  new  co-ordinates  are  orthogonal.  Equation  [4]  indi- 
cates that  for  fixed  t,  C and  6 differ  by  a constant,  hence  the  curves  of  con- 
stant C surfaces  of  constant  phase  while  the  curves  of  constemt  n point 

in  the  direction  of  r and  these  are  the  rays.  To  obtain  a unique  set  of  depen- 


dent variables  we  set 


K cos  $ , 


<2  * K sin  ^ 


where  ^ is  the  amgle  rays  make  with  the  x-axis.  Then  x and  y derivatives  in 
(81  and  [9]  are  converted  to  C and  n derivatives  with  ic,it'ta  as  dependent 


variables  . The  results  are 
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These  equations  (u:e  coupled  with  [10]  which  becomes 
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This  equation  is  used  to  eliminate  k from  [12] . 
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It  now  becomes  clear  how  the  term  geometrical  optics  applies.  In  linear 
optics  the  last  term  of  [13]  is  missing  and  if  the  second  derivatives  are 
negligible  as  well  then  ic»l.  Placing  this  into  [12]  we  have  iji^  = 0 so  that 
the  orientation  of  a ray  does  not  change  with  distance  along  the  ray.  Straight 
rays  lead  to  a geometrical  analysis.  The  term  geometrical  optics  approximation 
is  still  used  (see  II]  for  example)  to  refer  to  the  neglect  of  second  derivative 
terms  in  the  nonlinear  case  even  though  it  does  not  yield  straight  rays.  Also 
the  second  derivative  terms  are  still  referred  to  as  the  diffraction  terms  (as 
in  [1])  since  they  cause  rays  to  bend  even  though  the  nonlinear  term  also  causes 
bending  of  rays. 

Neglecting  the  diffraction  terms  of  [13]  we  eliminate  < from  [12]  obtaining 
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We  eure  able  to  notice  eui  exact  solution  by  interchanging  the  roles  of  dependent 
2md  independent  variables  by  a hodograph  transformation.  It  is  convenient  to  set 
p - a^  emd  then  we  have 

p _ . H-3p 

'p  ■ S'  S = ■ ■ 

One  can  now  eliminate  € or  n obtaining  a separable  equation  in  the  other.  A 
solution  that  represents  a focussing  beam  is 


[14] 


sin^ 


cosit 


A+2p  p/l+2p 

The  solution  written  relative  to  x and  y is  found  by  inverting  [11]  and  using 
[13]  (with  diffraction  neglected)  for  k . This  inversion  gives  x and  y as 
line  integrals  in  the  variables  t n then  [14]  converts  these  to  line 

integrals  in  ^ and  p . With  the  focal  point  placed  at  the  origin  to  fix  inte- 
gration constants.  Integration  gives 
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1 , 2p 

X = - ':r~  - log  ■; — — 
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cos  2$ 

2p(l+2p) 


[15] 


y = 


sin  2(t 
2p(l+2p) 


Equations  [15]  implicitly  determine  the  amplitude  P and  the  ray  orientation 
<(i  for  a suitable  domain  in  x and  y . We  establish  the  positions  of  rays  by 
fixing  n-  Setting  n = n in  the  first  equation  of  [14]  , iji  is  eliminated 
from  [15]  giving  rays  identified  by  n and  parameterized  by  p . The  resulting 
equations  are 
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h 


1+p 


p(l+2p) 


log 


l+2p 


i-(i+2p)n 

1+2  p 
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Similarly  one  can  obtain  a parameterized  form  of  the  surfaces  of  constant  phase 
by  fixing  6 = C in  the  second  equation  of  [14]  and  eliminating  0 from  [15] . 
The  rays  are  plotted  in  Figure  2a.  If  <ji  is  eliminated  from  [15]  one  can  obtain 
profiles  of  the  beam  at  various  positions  x and  these  are  plotted  in  Figure  2b. 


Solution  [15]  appears  in  (5)  and  further  focussing  solutions  appear  in  (6) 


A difficulty  with  these  solutions  is  that  in  the  region  of  the  focus,  equations 


[15]  fail  to  apply  since  the  expansion  for  small  a has  broken  down.  If  indeed 


there  is  a focus  where  the  amplitude  becomes  large  or  singular,  substitution  [3] 


is  unsuitable  because  the  oscillations  will  not  be  close  to  sinusoidal.  In  fact 


uniform  plane  wave  solutions  to  [1]  and  [2]  take  the  form  of  Jacobian  elliptic 


functions.  Clearly  [15]  cannot  be  correct  globally  but  since  the  only  indication 


of  trouble  is  in  the  focal  region  we  would  liope  that  [15]  holds  outside  the  focal 


region . Also  we  would  hope  that  it  is  the  geometrical  optics  approximation  that 


permits  the  singularity  to  occur  and  that  retaining  the  diffraction  terms  may 


prevent  the  focus.  This  is  investigated  in  the  next  section 


Aaplitude  Profiles  for  a Focussing  Solution 
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3 . Self -Trapped  Beams 

Having  noted  that  singularities  occur  in  the  geometrical  optics  approxima- 
tion we  return  to  equations  [B] , [9],  [10]  to  investigate  solutions  to  this 
fully  dispersive  version  of  the  problem.  Whenever  a local  maximum  of  amplitude 
occurs,  nonlinearity  without  diffraction  causes  focussing.  If  we  are  now  able 
to  produce  solutions  where  rays  converge  initially  due  to  a locally  large  ampli- 
tude, but  do  not  further  on  produce  a focus,  we  shall  have  strong  evidence  that 
diffraction  is  an  essential  effect  in  nonlinear  optics. 

We  first  find  a solution  that  represents  a uniform  self-trapped  beam  by 
seeking  solutions  independent  of  x so  that  rays  are  straight  and  nonlinearity 
and  diffraction  are  balanced.  Taking  x as  the  direction  of  propagation,  if 
the  properties  of  the  beam  do  not  vary  in  x then  is  constant,  is 

zero  and  a depends  only  on  y . Already  [8]  and  19]  are  satisfied  and  setting 
= K > Kq,  a solution  to  110]  is 


A^-1 


/-2 

sech  / K -1  y 


This  solution  decays  for  large  [y]  while  all  other  solutions  oscillate  in  y . 
This  solution  has  been  discussed  by  Chiao  (4)  , AWmanov  (1)  and  many  others  al- 
though most  of  them  find  it  as  a solution  to  approximate  equations  since  it  satis- 
fies their  paraxial  restrictions. 

To  upset  the  balance  that  maintains  straight  rays  we  perturb  this  solution 
by  functions  that  vary  slowly  in  x . We  set  x - ex,  where  e is  a small  posi- 
tive constant,  and  seek  solutions  to  [8] , 19) , IlO)  that  are  power  series  in  e 
of  the  form 


117]  . 


a - a(y)  + cGj^(x,y)  + ... 

= ic  + e0j^(x,y)  + ... 

Kj  - eVj^(x,y)  + e^V^(x,y)  + ... 


As  it  happens  this  expansion  gives  acceptable  solutions  for  the  quantities 

displayed  but  to  extend  to  higher  accuracy  the  perturbations  have  a small  effect 

on  the  original  solution  a(y}.  This  is  accounted  for  by  replacing  y by  the 

2 3 

new  variable  y=  (1+  eYj+  •••)y  where  the  are  determined  by  sup- 
pressing secular  terms  in  the  at  each  stage  of  solution.  For  simplicity  we 

shall  only  seek  beams  that  are  symmetric  about  the  x-axis.  This  requires  that 
and  be  even  in  y and  that  and  be  odd-  in  y . 

Placing  117]  into  [8] , [9] , [10]  and  equating  coefficients  of  like  powers 
of  e we  obtain  a hierarchy  of  equations.  From  [8]  we  obtain  «ui  equation  that 
determines  to  be 

= C cosh^  / ic^-l  y 
and  the  equation  determining  is 

[18]  (aV)  + {2aG-V  1 + 2aKG  - + a^U  - = 0 . 

2 y 1 1 y lx  lx 

From  [9]  we  find  that  must  be  a function  of  x , hence  we  set 

= u(x)  . 

Finally  [10]  gives  the  equation  determining  Gj^  as 

[19]  Gj^yy  - (k^-I)Gj^  - 2x30 + 6a^Gj^  = 0 . 

Since  must  be  odd  we  take  C = 0 and  this  necessitates  finding 

to  obtain  the  y-component  of  x . While  is  an  arbitrary  function  of  x , 

succeeding  terms  in  expansion  [17]  will  not  be  small  unless  u(x)  and  all  its 
derivatives  are  slowly  varying.  Hence  u(x)  must  be  analytic  in  x.  When 
is  placed  into  [19]  the  x-dependence  factors  out.  Setting  Gj^  = u(x)g(y),  [19] 
becomes 

-2  2 / —2  - / -2  /_2 
- (k  -1)  (1-6  sech  V K ~X  y)g  = 2k/  k -1  sech  / k -1  y - 

The  bpunded  solution  of  this  equation  is  found  by  differentiating  [16]  with 

respect  to  k and  y since  perturbing  k by  a constant  gives  another  solution 

of  the  form  [16] . Thus  we  obtain 


1 


I 


I 
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(20] 


I 


= u(x) 


c^-l 


sech  /k^-1  y - k y sech  /k^  -1  y tanh  /ic^  -1 


+ A sech 


y tanh 


y J 


Since  Gj^  must  be  even  in  y we  take  A » 0 . Then  we  solve  118]  for 
obtaining 


[211 


K^-1 


duci)  r -2 

dx 


r y + 


2k^-1 


2(x^-l)^'^^ 


slnh  2 / K^-1  y 


The  perturbation  is  now  fixed  when  the  arbitrary  function  u(x)  is  specified. 

To  obtain  the  beam  of  Figure  lb  we  could  take  u(x)  » cos  x , for  example, 
to  have  a beam  with  periodically  converging  and  diverging  rays  but  no  focus. 

we  could  take  u(x)  * x to  have  a beam  that  does  focus.  Thus 
we  assert  that  a focus  may  or  may  not  result  when  waves  in  a nonlinear  focussing 
medium  begin  to  converge  and  that  diffraction  is  sufficient  to  prevent  focussing 
although  it  does  not  necessarily  do  so. 
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